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Abstract. We used the photometric database of the Sloan Digital Sky Survey (SDSS) to estimate the reddening 

of 13 SDSS quasars selected on the basis of the presence of zinc absorption lines in an intervening Damped Ly a 
ivj ', (DLA) system. In 5 of these quasars the reddening is detected at > 2a confidence level in two independent color 

^+ ' indices of the SDSS ugriz photometric system. A detailed analysis of the data supports an origin of the reddening 

f^ ' in the intervening absorbers. We used these rare measurements of extinction in DLA systems to probe the relation 

^O . between extinction and metal column density in the interval of absorption redshift 0.7 < z < 2.0. We find that 

^^ ' the mean extinction in the V band per atom of iron in the dust is remarkably similar to that found in interstellar 

fH ^ clouds of the Milky Way. This result lends support to previous estimates of the dust obscuration effect in DLA 

Oj' systems based on a Milky Way extinction/metal column density relation. We propose a simple mechanism, based 

on dust grain destruction/accretion properties, which may explain the approximate constancy of the extinction 

per atom of iron in the dust. 
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S^ 1. Introduction From Milky Way interstellar studies it is well known 

H that the extinction per H atom is roughly constant, with 

Ci , Interstellar dust is a pervasive component of galaxies and ^ ^^^.^^^ ^^^^^ ^^/^^ ^ 5 3 ^ ^o_22 ^^g ^^2 ^j^^^^-^ 

plays a key role in a variety or astrophysical processes ai- , , in-ToN u 4 • au i- j-- • j-u ta u j 

'^ ■^ ■' ■' r 1 et al. 1978), where Ay is the extmction m the V band, 

lecting, with its presence, many types or observations or rpi i i- • i-j -j • j. i c i j 

, ' , , , , , r /r. ^'^^ relation is valid over a wide interval oi column den- 

the local and high-redshift Universe (Spitzer 1978, Drame ^.^.^^^ .^^^ ^ ^ ^^ig < ^^[^^-2^ < 3 ^ iq^i. Observations 

2003, Meurer 2004). Understanding whether or not the ^ , , , ?" /„ , r , '^ 1 r.^^^\ ■ ,• 

, ^ ^ , . , , , Ti 1 ■ -1 i- ■^^ ^^ ^ 01 translucent clouds (Kachtord et al. 2002) indicate that 

dust at high redshiit shares similar properties with that ,, , . , , , ; , . , , „ ,. ,. , 

, , 1 1 . r. . ; 1 ; , 1 . ; ; . .1 1 the sauic rclatiou holds at high values oi extinction up to 

at low redshitt IS central to correctly interpreting the ob- . „ \ i, . • , , i r ^ /^r ■ i 

r , . , . ,, 1 TT • At/ ~ d mag , albeit with a modest trend or Ai//A'h with 

servations ot galaxies and quasars m the early Universe. □ _ /i // ^ a \ rr^, ■ nnno\ 

In this paper we focus our attention on a particular 



Rv = Av/{Ab - Av) (Draine 2003). 



property of the dust, namely the extinction per unit col- ^he fact that Ay scales linearly with N^ over a wide 

umn density of metals. The aim of this work is to esti- "^ange of physical conditions probably reflects an underly- 

mate this quantity in individual clouds at high-redshift ^^S property of dust grains. If this property is universal, 

and compare it with measurements performed in local ^^ expect to find its signature m high redshift clouds, 

clouds. The results of this investigation are relevant for To investigate this possibility we must understand which 

studies of the extinction generated by quasar absorption ^^ ^^^ physical relation underlying the empirical relation 

line systems and for casting light on the physics of dust ^v oc N^. Since the extinction is generated by the dust 

grains at different cosmic epochs. particles. Ay should scale with the column density of 

atoms in the dust rather than with the column density 

Send offprint requests to: G. Vladilo of hydrogen in the gas. For Galactic clouds, it is easy to 
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show that Nn oc N-n, where N-n is the dust-phase column 
density of refractory elements^. Therefore, the empirical 
Milky Way relation Ay ex iVn might be the result of a 
more universal relation Ay ex TVr. The scientific goal of 
our investigation is to test the existence of this relation 
from a study of high redshift interstellar clouds. If the re- 
lation exists, it should also be able to predict the values 
of Ay /N-ti observed in the Magellanic Clouds, which are 
lower than in the Milky Way, even if relatively constant 
within each galaxy (e.g. Gordon et al. 2003). 

To accomplish our goal, wc need high redshift mea- 
surements of extinction and column densities, which are 
rather difficult to obtain. To measure the column densities 
one needs a bright, point-like source located beyond the 
cloud. Among sources at cosmological distances, QSOs are 
the best targets given the very fast decline of Gamma Ray 
Bursts. To measure the extinction one needs to estimate 
which fraction of QSO light has been removed from the 
line of sight due to absorption and scattering processes 
within the cloud. This task is particularly difficult owing 
to our poor knowledge of the true spectral energy distri- 
bution of individual QSOs which, in addition, is affected 
by variability. 

A way to overcome these difficulties is to compare com- 
posite spectra of quasars with and without foreground ab- 
sorption systems in a search for a systematic change of the 
continuum slope. After the pioneering work by Pei et al. 
(1991) and the lack of reddening signal from DLA systems 
at redshift z ~ 2.8 (Murphy & Liske 2004), this approach 
led to the detection of reddening in front of Ca ii systems 
(Wild & Hewett 2005, Wild et al. 2006) and Mgii systems 
(Khare et al. 2005, York et al. 2006) at lower redshift. 

In the present work we do not follow the same ap- 
proach, since we are interested in individual absorbers 
rather than in statistical samples. In particular, we are 
interested in studying DLA systems, the QSO absorbers 
with highest Hi column density, A^(IIi) > 2 x 10^° atoms 
cm~^, that trace the interstellar medium of high redshift 
galaxies (Wolfe et al. 2005). Our goal is to estimate the 
reddening of individual DLA systems by means of photo- 
metric techniques that can be efficiently implemented in 
large databases. 

Studies of photometric reddenings of individual 
quasars with foreground DLA systems have started only 
recently (Khare et al. 2004, Ellison et al. 2005). Two of 
Ellison et al.'s (2005) DLA quasars seem to be significantly 
reddened, but in one of these two cases the reddening may 
well be due to dust in the quasar host galaxy. A prelim- 
inary study of the reddening versus gas-phase depletion 



^ The symbol A'' is introduced to indicate column densities of 
atoms in the dust phase. The validity of the relation A''h oc A'^r, 
is the result of the constant (solar) composition of the local 
interstellar medium. If the composition is constant, the gas- 
phase column density of volatile elements scales linearly with 
the dust-phase column density of refractory elements (Vladilo 
& Peroux 2005). Hydrogen is mostly present in the gas phase 
and is, in practice, a volatile element. 



of the absorber, based on a very small sample, has been 
reported by Khare et al. (2004). The behaviour of the 
rest frame extinction versus metal column density of the 
absorber, the object of the present study, has not been 
investigated. 

In the previous photometric studies the intrinsic color 
of the quasars was estimated using a composite spectrum 
of quasars obtained from the Sloan Digital Sky Survey 
(SDSS) database (e.g. Richards et al. 2001). The error as- 
sociated with the composite spectrum was not propagated 
into the reddening measurement. 

In the present work we measure the reddening us- 
ing the quasar photometric database of the SDSS (e.g. 
Schneider et al. 2005), paying special attention to the 
treatment of the errors. Instead of using a quasar compos- 
ite spectrum, we build up a distribution of quasar colors 
at each redshift to estimate both the intrinsic color (the 
median of the distribution) and its error (the dispersion of 
the distribution). This error is then propagated into the 
reddening measurement. At variance with previous work, 
we estimate the reddening of each quasar in at least two 
different color indices. Using different colors yields two ad- 
vantages. First, we can minimize local systematic errors 
(e.g. contamination of the quasar continuum in one spe- 
cific passband of the spectrum). Second, as wc show below, 
we can test if the reddening originates in the intervening 
system. 

In the selection of the targets we only consider quasars 
with a single DLA system. In this way we can unambigu- 
ously investigate the relation between quasar extinction 
and column densities of the absorber. In the next section 
we describe the sample, while in Sect. 3 we present the red- 
dening measurements and their conversion to rest-frame 
extinction. In Sec. 4 the extinction is compared to the 
dust-phase column density of the refractory element iron, 
A^Fo, to test the existence of a general relation Ay oc A^r 
in high redshift and local interstellar clouds. The results 
are discussed in Sect. 5 and the whole work is summarized 
in Sect. 6. 

2. The sample 

In order to generate our sample we followed three crite- 
ria. The first criterion was driven by the requirement of 
estimating the dust-phase column density of iron. As we 
shall see below, this can be obtained from a simultaneous 
measurement of the zinc and iron column densities. We 
therefore gathered all the quasar absorption systems with 
detections of Znii and Fen lines. 

Only absorption systems originating in neutral regions 
with A(Hi) > 2 X 10^° atoms cm~^, i.e. damped Lyman 
a (DLA) systems, were considered. In these regions, Znii 
and Fe ii are expected to be the dominant ionization stages 
of zinc and iron and the column densities of these metals 
can be determined without applying ionization corrections 
(see Vladilo et al. 2001 and refs. therein). 

In the selection process wc used our previous compila- 
tions of DLA systems (Vladilo 2004, Kulkarni et al. 2005) 
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supplemented by results from recent work (Wang et al. 
2004, Akerman et al. 2005, Peroux et al. 2006, Meiring et 
al. 2006), collecting a total of about 70 absorbers. 

The second criterion was to choose, from the resulting 
list, only the DLA systems with the quasar included in the 
SDSS photometric database. We used in most cases the 
Data Release 3 (DR3) catalog by Schneider et al. (2005), 
supplemented by the Data Release 4 (Adelman-McCarthy 
et al. 2006) when necessary. For a bright QSO not included 
in these catalogs (see Table 1) we used the list published 
by Richards et al. (2001). 

The homogeneous database of SDSS photometry is an 
ideal tool for building up a control sample of quasar colors 
to be used in the reddening measurement. Thanks to the 
large size of the SDSS catalog one can easily obtain a con- 
trol sample with large numbers even in very narrow bins 
of redshift and magnitude. The presence of simultaneous 
photometry in the 5 ugriz filters allows us to obtain re- 
liable reddening measurements in different color indices. 
Estimates of the mean absorber reddening based on pho- 
tometric measurements have been shown to be reliable by 
comparison with those based on spectroscopic measure- 
ments (York et al. 2006). 

The third criterion was to exclude lines of sight with 
multiple DLA systems. In fact, for these cases the associa- 
tion between quasar reddening and metal column density 
would be ambiguous. Lines of sight with one neutral re- 
gion and one (or more) regions of high ionization (e.g. 
C IV, Siiv) were kept in the sample, assuming the redden- 
ing of the ionized region(s) to be negligible compared to 
that of the neutral region. 

As a result of this selection, we obtained the sample of 
13 pairs of QSOs/DLA systems shown^ in Table 1. The 
large rate of rejection relative to the original Znii sample 
is probably due to several factors. One could be the poor 
overlap in magnitude space between the quasars of the 
Znii sample, mostly bright (m < 18), and those of the 
SDSS catalog, mostly fainter. Another could be the limited 
sky coverage of the SDSS. Only a very few SDSS quasars 
with detected Zn were rejected because of the presence of 
multiple DLA systems. 

The basic data for the selected absorption systems are 
given in the last 5 columns of Table 1. The original spec- 
troscopic studies of the metal lines were based on observa- 
tions collected with the ESO VLT telescope (4 cases), the 
Multiple Mirror Telescope (4 cases), the Keck telescope 
(3 cases), the 4-m telescope at the Kitt Peak National 
Observatory (1 case) and the HST (1 case). 

For most of the selected systems the DLA nature is 
confirmed by a direct measurement of the damped Lya 
profile. At redshift 2;abs > 1-9 the A^(Hi) measurements 
are based on observations obtained at the ground-based 
facilities mentioned above. The Lya measurements at 
Zahs < 1.9 are based on UV observations (HST in 6 cases 
and lUE in 1 case). For two systems at Zabs < 1-9 without 



UV observations, the DLA nature is suggested either by 
the unusually strong Znii lines (J0121-I-0027) or by the 
presence of many species of low- ionization, including Mg i 
(J2340-0053). 

As a final check of the sample of Table 1 , we searched 
for absorption lines indicative of potential sources of red- 
dening, even if not classified as DLA systems. In this 
search we used the low-resolution SDSS spectra of the 
quasars, plus all the information presented in previous 
studies of high resolution spectra. A direct search for 
strong Lya lines at z < Zqso was not particularly use- 
ful since for all except one of our targets the portion of 
Lya forest covered by the SDSS spectrum is too narrow 
or non-existent. However, the wavelength coverage is suf- 
ficiently large for detecting Mg ii absorptions over a large 
redshift interval (0.36 < Zabs < 2.29). We expect that any 
intervening H i cloud with high extinction would produce 
a strong Mgii resonance doublet. In fact, it is now proven 
observationally that Mg ii absorbers do contribute to the 
reddening of the quasars (York et al. 2006) and therefore 
can be used to trace additional sources of dust. We there- 
fore searched for narrow, but intense absorption features 
that could be attributed to Mg ii systems at a redshift dif- 
ferent from that of the selected DLA system. In this search 
we also checked for Mgii absorptions superposed on the 
Mg II quasar emissions, an indirect signature of dust in the 
quasar host galaxy. In most cases we found no signature of 
strong Mgii lines or other strong lines from Hi regions, ei- 
ther from intervening systems or from associated systems 
(quasar host galaxy). 

In two quasar spectra, however, we did find evidence 
for additional absorption from neutral gas. These quasars 
are listed separately at the end of Table 1 since for these 
two lines of sight an additional contribution to the redden- 
ing might be present, in addition to that expected from 
the DLA system. Details on the additional absorptions are 
reported in the footnotes to the table. With the possible 
exception of these two cases, the DLA systems of Table 
1 are very likely to be the major source of reddening of 
their respective quasars. 

Given the pre-selection in Znii and Fen, the sample 
of Table 1 does not represent a random sub-set of the to- 
tal DLA population. Given the current limitations of high 
resolution spectroscopy at A > 800 nm, the pre-selection in 
Znii precludes the range of absorption redshift z > 3. In 
spite of these limitations, the list of Table 1 is the best cur- 
rently available sample for studying possible correlations 
between reddening and metal column densities of individ- 
ual systems, considering that reliable Znii and Fell col- 
umn densities are extremely difficult to derive from SDSS 
data for individual objects. 



The SDSS identifier given in Table 1 is truncated to 4 digits 
in RA and DEC in the rest of the paper. 
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Table 1. SDSS quasars with previously detected Znii and Fell absorptions in an intervening Damped Lya system. 



SDSS 


Other name 


■Sqso 


9 

(mag) 


■^abts 


logAr(Hi) 
(cm-2) 


Ref." 


[Zn/H]" 


Ref.<= 


J001306. 1+000431 


LBQS 0010-0012 


2.165 


18.65 


2.025 


20.80 ±0.10 


LPS03 


-1.05 ±0.10 


LPS03 


J001602.4-001225 


Q0013-004 


2.087 


18.28 


1.973 


20.83 ±0.05 


PSL02 


-0.65 ± 0.06 


PSL02 


J012147.7+00271 


B0119+0011 


2.224 


20.00 


1.388 


[21.04]'' 


W04 


— 


W04, P05 


J093857.0+412821 


Q0935+417 


1.936 


16.49 


1.373 


20.52 ±0.10 


LWT95 


-0.92 ±0.11 


MLW95 


J094835. 9+432302 


Q0948+433 


1.892 


18.10 


1.233 


21.62 ±0.06 


R05 


-1.12 ±0.06 


R05 


J101018. 1+000351 




1.400 


18.27 


1.265 


21.52 ±0.07 


R06 


-1.02 ±0.09 


MOO 


J110729. 0+004811 




1.391 


17.64 


0.741 


20.98 ±0.15 


R06 


-0.60 ±0.16 


K04 


J115944. 8+011206 


Q1157+014 


2.000 


17.59 


1.944 


21.80 ±0.10 


WB81 


-1.36 ±0.13 


PSLOO 


J123200. 0-022404 


PKS 1229-021 


1.044 


17.13 


0.395 


20.75 ± 0.07 


B98 


-0.47 ±0.14 


B98 


J132323. 7-002155 = 




1.388 


18.46 


0.716 


20.21 ±0.20 


P06, R06 


+0.57 + 0.21 


K04, P06 


J150123. 4+001939 




1.928 


18.11 


1.483 


20.85 ± 0.05 


R06 


-0.40 + 0.07 


MOO 


J223408.9+000001-f 


LBQS 2231-0015 


3.015 


17.57 


2.066 


20.56 ±0.10 


LW94 


-0.75 ±0.10 


PW99 


J234023. 6-005326" 




2.085 


17.76 


1.361 


[20.30]'' 


K04 


— 


K04 



" References for H I column density data. 

We adopt the usual definition [X/Y] = \og{Ny^/Ny-) — log(X/Y)Q. Throughout this paper we use the meteoritic solar abundances of Anders Sz 
Grevesse (1989) for consistency with most previous work on interstellar depletions. 
^ References for Zn II and Fell column density data (see also Table 3). 

Indirect estimate of log A'^(Hl) published by the authors based on reddening determinations. 
^ For the absorber at ^abs — 0.716 towards this quasar the Hi column density adopted here (Pcroux et al. 2006) is lower than that given by Rao et 



al. (2006), logA^(Hi) ^ 20. 54^ 



The adopted value formally lies below the DLA definition threshold, but we keep this system in the list because 



the diff'erence (ci 0.1 dex) is significant only at c=i 0.5 cr level and because, from a study of the ionization properties of this absorber (Peroux et al. 

2006), we do not find differences relative to the properties typical of DLA systems. 

■' The SDSS photometric data for this bright quasar (Foltz et al. 1989) can be found in Richards et al. (2001). No SDSS spectrum is available to 

date. A large number of absorption features is present in the spectrum published by Lu & Wolfe (1994); in addition of the DLA system at z — 2.066, 

these authors report the presence of a strong system with both neutral and ionized gas at z — 2.6527. 

^ The SDSS spectrum shows a large number of absorption features, including a Mg II doublet at z^hs — 2.05; in spectra of higher resolution Khare 

et al. (2004) find a system at 2 — 2.0547 with a mix of ionization states, including C I, a signature of neutral gas. 

References. B98: Boisse et al. (1998); K04: Khare et al. (2004); LW94: Lu & Wolfe (1994); LWT95: Lanzctta et al. (1995); LPS03: Ledoux et al. 

(2003); M06: Meiring et al. 2006; MLW95: Meyer et al. (1995); PSLOO: Petitjean et al. (2002); PSL02: Pctitjcan et al. (2002); P06: Pcroux et al. 

(2006); P05: Prochaska (2005, priv. comm.); PW99: Prochaska & Wolfe (1999); R05: Rao et al. (2005); R06: Rao et al. (2006); W04: Wang et al. 

(2004); WB81: Wolfe & Briggs (1981). 



3. Reddening measurements 

For each quasar of Tabic 1 we measured the reddening in 
the color index (y — x) from the expression 



My - ^) = (y -x) - (y - ^)c 



(1) 



where (y — x) is the observed color of the quasar and 
{y — x)o the intrinsic color in absence of reddening. All 
quantities in this definition are in the observer's frame. 

The color {y — x) was measured using different pairs 
of PSF magnitudes x,y = u,g,r,i,z of the SDSS DR3 
catalog preliminarily corrected for Galactic extinction 
(Schneider et al. 2005). 

We only used bandpasses falling in spectral regions 
where the continuum of the quasar is not contaminated by 
the Ly a forest or by the quasar Lya emission. In practice, 
taking into account the effective wavelengths and FWHM 
of the ugriz bandpasses (Fukugita et al. 1996, Schneider 
et al. 2005), this criterion precludes the use of the u band 
when Zqso > 1.67 and also of the g band when Zqso > 2.25. 
As an additional criterion, pairs of adjacent ugriz band- 
passes were not considered since wavelength separation 
is essential for detecting the reddening of the quasar, if 
present. The colors obtained with these criteria are given 
in Table 2, Col. 3. 

To estimate the intrinsic color and its uncertainty we 
first generated a control sample of reference quasars for 
each quasar of our list. Homogeneity in redshift and ap- 
parent magnitude were the criteria adopted for building 
each control sample. 



The homogeneity in redshift is very important given 
the redshift dependence of the quasar colors (Richards et 
al. 2001). The large database of the DR3 catalog (46420 
objects) allowed us to select control samples with a large 
number of reference quasars even using small redshift bins. 
In practice, we adopted z^ci G (zqso — 0.05, Zqso + 0.05) in 
all cases but one (see note a of Table 2). 

After selecting a sample with a given redshift, we ap- 
plied the criterion of homogeneity in apparent magnitude. 
This is in practice equivalent to a homogeneity in abso- 
lute luminosity since the quasars of each sample are at the 
same redshift. We used the infrared magnitude z, which 
is the least affected by extinction. The width of the mag- 
nitude bin was tuned in such a way as to obtain control 
samples of the same size for all the quasars of our list. 
For instance, to obtain a control sample of 100 quasars, 
starting from the sample already binned in redshift, the 
typical width of the magnitude bin was of « 2 magnitudes. 
Changes of the intrinsic slopes of the quasar continua over 
an interval of « 2 in absolute magnitude are expected to 
be very modest (e.g. Yip et al. 2004). 



For each control sample and color index of interest 
we then derived the distribution of colors corrected for 
Galactic extinction. The median of the color distribution 
was adopted as an estimator of the intrinsic color (y — x)o. 
Compared to the mean, the median offers the advantage 
of being very stable for different choices of the control 
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Table 2. Colors and reddenings of the quasars of Table 1. 



SDSS 


Color 


Quasar color 




Median 




PR" 


Quasar reddening 


^qso 


index 


(v-^) 


3300° 


(y - x)o 
sioo' 


Clean^ 


(%) 


A(y-x) 


J1232-0224 


(u-z 


0.515 ± 0.023 


0.323 


0.324 


0.325 


6 


+o.i9ii;il 


1.044 


{u-i^ 


0.551 ± 0.020 


0.322 


0.322 


0.323 




+0.23t«„-= 




(g~z 


0.238 ± 0.029 


0.173 


0.177 


0.178 




+0.06tl\l 




(«-r 


0.572 ± 0.020 


0.365 


0.362 


0.361 




+0.2ll«:- 




(g-i) 


0.275 ± 0.027 


0.176 


0.180 


0.180 




+o.o9iS:ii 




(r~z 


-0.058 ±0.020 


-0.038 


-0.016 


-0.018 




-o.o4i«:- 


J1323-0021 


[u — z 


1.322 ±0.033 


0.361 


0.394 


0.387 


6 


+0-9it'oil 


1.388 


{u — i 


1.249 ± 0.033 


0.385 


0.400 


0.390 




+0.86tlil 




(s-^ 


0.827 ± 0.033 


0.255 


0.289 


0.293 




+0.531-.- 




{u — r 


1.037 ± 0.028 


0.315 


0.325 


0.320 




+0 72+"-- 




(s-i) 


0.754 ± 0.033 


0.288 


0.294 


0.293 




+0A6t"„\l 




{r-z 


0.285 ±0.026 


0.033 


0.059 


0.055 




+o-^3tl:ll 


.11107+0048 


(u~z 


0.459 ± 0.029 


0.351 


0.364 


0.335 


11 


+ 0.12_o 2g 


1.391 


{u~i] 


0.495 ± 0.024 


0.383 


0.401 


0.377 




+0 12+"-^'^ 




{g-z 


0.342 ± 0.027 


0.256 


0.261 


0.238 




+o.iol"„:S 




{U-T 


0.413 ±0.026 


0.314 


0.324 


0.305 




+0.11;-- 




(g~i) 


0.378 ± 0.022 


0.290 


0.281 


0.270 




+0 11+011 




{r~z 


0.046 ± 0.026 


0.037 


0.026 


0.018 




+o.o3+;;f, 


JlOlO+0003 


[u — z 


0.500 ± 0.037 


0.337 


0.388 


0.374 


5 


+0.13+-- 


1.400 


(u^V 


0.520 ± 0.027 


0.356 


0.393 


0.380 




+0.141S:- 




(g-z 


0.159 ±0.032 


0.244 


0.253 


0.250 




-omtlil 




(u-r 


0.486 ± 0.025 


0.294 


0.315 


0.311 




+0.l8tlit 




(g-i) 


0.180 ±0.019 


0.279 


0.289 


0.283 




-o.iots:!! 




{r-z 


0.014 ±0.032 


0.039 


0.039 


0.032 




-0.02i-.f3 


J0948+4323 


(g-z 


0.289 ± 0.057 


0.365 


0.374 


0.339 


29 


-0.051-.- 


1.892 


(g-i) 


0.397 ± 0.054 


0.303 


0.302 


0.278 
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" Median eolor of tiie eontrol sample of tiie 300 quasars with closest value of z magnitude. 
Half width of redshift bin Sz = 0.05 except for J2234+0000 (Sz = 0.15). 

Median color of the control sample of the 100 quasars with closest value of z magnitude. 
Half width of redshift bin Sz — 0.05 in all cases. 
^ Median color of the 'Clean' sub-sample of SlOO (the quasars without absorption features in their SDSS spectra). 

Percent of quasars rejected from the control sample SlOO on the basis of the presence of absorption lines in their spectra. 



sample, even when some quasars with anomalous colors 
happen to be included in the sample. 

We performed several tests to assess the robustness of 
the median to sampling errors. The results of one of these 
tests is shown in Cols. 4 and 5 of Table 2, where we com- 
pare the medians of the control samples labelled 'S300' 



and 'SlOO', with 300 and 100 reference quasars respec- 
tively. One can see that the differences in the medians are 
in most cases < 0.01 magnitudes. Other tests performed 
using sub-samples of S300 and SlOO with statistically sig- 
nificant number of objects yield the same indication for 
the magnitude of the sampling error. 
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Since a fraction of the quasars of the control sample 
might be affected by reddening, we made a visual inspec- 
tion of the SDSS quasar spectra in the 'SlOO' samples 
to search for signatures of intervening absorption systems 
which could, in principle, redden the reference quasars. 
In most quasars of our list we could not search for inter- 
vening Damped Lyman a absorptions since the Lyman 
a forest lies out of the observable wavelength range. We 
therefore searched for strong, narrow absorption lines red- 
wards of the forest, tentatively identified as species of low 
ionization that might arise in an intervening Hi region. 
Quasars showing absorption lines with residual intensity 
< 50% redwards of the quasar Lyman a emission were 
rejected from each SlOO sample. Also Broad Absorption 
Line (BAL) quasars were rejected. In this way we gen- 
erated a new control sample, labelled 'Clean' in Table 2. 
The fact that the percentage of rejection tends to increase 
with the quasar redshift (see column labelled 'PR') argues 
in favour of an origin of the rejected absorptions in inter- 
vening regions. When the percentage of rejection is sig- 
nificant, one can see that the median color of the parent 
sample SlOO is systematically redder than the correspond- 
ing median of the 'Clean' sub-sample. Even if the effect is 
modest, and often of the same order of the sampling er- 
rors, this result indicates that the intervening absorptions 
do contribute to the reddening of the quasars. Also the 
dispersion of the color distribution shows a little, but sys- 
tematic change, being smaller after the rejection process. 
We therefore adopted the color distribution of the 'Clean' 
control sample for our measurements. 

The intrinsic color [y — x)o and its error were estimated 
from the median and the dispersion of the 'Clean' distri- 
bution, respectively. Even though the dispersion around 
the median is commonly estimated using the interquar- 
tile range, i.e. the interval between the 25th and the 75th 
percentiles, we preferred to adopt the range between the 
16th and the 84th percentiles, which is more conservative 
and consistent with the range bracketed by ±1 standard 
deviations in a normal distribution. In fact, by definition, 
the adopted interval brackets 68% of the area of the dis- 
tribution around its central value, in the same way as the 
range ilcr does for a normal distribution. 

The resulting reddening is given in the last column 
of Table 2 for each quasar and color index. The adopted 
uncertainty of the reddening measurement was obtained 
by propagating the error of the quasar color (column 3 
of the table) and the dispersion of the color distribution 
explained above. One can see from Table 2 that quasar 
reddening is detected at the ~ 3ct level in two objects and 
at the ~ 2(7 level in other three objects. 

A test of the accuracy of the zero point of our measure- 
ments is offered by the quasar J0013-I-0004, whose fore- 
ground DLA system is dust-free on the basis of its solar 
[Zn/Fe] ratio (see Table 4 and Section 4). For this quasar 
we find zero reddening within the errors in all the color 
indices considered, as expected for a dust-free system. 
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Fig. 1. Comparison of different types of extinction curves 
normalized to the V band, ^(A) = A\/Av- MW curves: 
CCM model adopted in this work (thick line) versus P92 
model (dotted line). SMC curves: modified P92 model 
adopted here (thin line) to fit the G03 data (circles) versus 
original P92 model (dashed line). See Section 3.1. 

3.1. Conversion to rest-frame extinction 

The quasar reddening measured in the observer's frame, 
A(y — x), can be converted to the V band extinction in 
the rest- frame of the absorber. Ay, if one knows the nor- 
malized extinction curve of the absorber, ^(A) = Ax/ Ay, 
and its redshift, Zabs. In the Appendix A we show how to 
perform this conversion taking into account the slope of 
the quasar continuum. If the bandpasses are sufficiently 
narrow, the variation of the continuum accross the band- 
pass can be neglected. In this case, we obtain the simple 
relation useful for the purpose of discussion. 



Av 



A(2/ 



a 



X-u 



A a; 



14-2ab. 



■)-C(l^) 



(2) 



where A^;^ and Xy^ are the wavelengths of maximum re- 
sponse of the bandpasses x and y. The values of Ay pre- 
sented in our work, however, were derived from the more 
general relation l|A.10(l which is also valid for wide band- 
passes. The response curves of the SDSS ugriz bandpasses 
were taken from Fukugita et al. (1996 and priv. comm.). 
We adopted a power law with spectral index a^ = —0.5 
for the quasar continuum. 

In our computations we used the two extinction curves 
most commonly adopted in the literature, i.e. the average 
curve of the Milky Way, characterized by an extinction 
bump at 217.5 nm and that of the Small Magellanic Cloud 
(SMC), characterized by a fast UV rise and without the 
bump (see e.g. Draine 2003). For the first one we adopted 
the model of Cardelh et al. (1988; CCM) with Ry = 3.1. 
For the SMC curve we fitted the model of Pei (1992; here- 
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Fig. 2. Extinction Av of the reddened quasars of Table 2 estimated in the rest frame of the intervening DLA system 
and plotted versus color index for different extinction curves. Color index coding: 1 = (u — z), 2 — (u^i), 3 = (g — z), 
4 = (m — r), 5 = (g — i); 6 = (r — z). Circles: SMC extinction curve. Diamonds: MW extinction curve. See Section 3.1. 



after P92) to the average SMC bar data of Gordon et al. 
(2003; hereafter G03). As we show in Fig. ^ the adopted 
MW and SMC curves are very similar to the models of 
P92 often used in the literature. 



From previous studies we expect the incidence of MW 
curves to be very low, compared to that of SMC curves, 
among quasar absorbers (Wild & Hewett 2005; York et al. 
2006 and refs. therein). However, in a few cases evidence 
for MW-type curves has been reported, with one clear de- 
tection of the bump in a DLA system at Zabs — 0.524 ob- 
served with the HST (Junkkarinen et al. 2004), and three 
possible identifications in absorption systems at Zabs — 1-4 
- 1.5 observed in SDSS quasar spectra (Wang et al. 2004). 



In principle any type of extinction curve can be used 
in the conversion from reddening to rest-frame extinction. 
In practice, however, the results obtained from flat ex- 
tinction curves, such as the MW curves with Ry > 5, are 
unreliable. This can be understood from the simple rela- 
tion |2Il: if the extinction curve is exactly flat the measured 
A[y — x) is divided by zero in the conversion. 



3.2. Comparison of results obtained from different 
color indices 

The rest frame extinction Ay is an intrinsic property 
of the absorber and must be independent of the band- 
passes X and y used in the measurement of the redden- 
ing. Therefore, by comparing values of Ay obtained from 
different measurements of A(j/ — x) in the same quasar 
we can test the hypothesis that the measured reddening 
originates in the intervening system. If the hypothesis is 
correct, we must obtain a constant Ay from the different 
values of A(y — x) converted using the redshift Zabs and 
the extinction curve ^(A) of the absorber. 

The result of this test depends on the adopted extinc- 
tion curve. However, all known extinction curves share a 
very similar slope in the spectral region redwards of the 
217.5 nm bump (see Fig. ^. Therefore, when z^hs is suffi- 
ciently low for all the bandpasses to fall redwards of the 
redshifted bump the test is, in practice, independent of 
the adopted curve. 

In Fig. [21 we plot Ay versus color index for the av- 
erage MW and SMC extinction curves. Only quasars for 
which the reddening was detected in at least three colors 
are shown in the figure. One can see in Fig.[21that in each 
quasar Ay is constant within the errors for at least one of 
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the two extinction curves. The test is particularly strin- 
gent for J1323— 0021, for which we have at our disposal 6 
different measurements. These results are consistent with 
the hypothesis that the reddening originates in the inter- 
vening DLA systems. 

As expected, the test is unable to discriminate between 
different types of extinction curves for the absorbers at 
lower redshift. However, when Zabs is sufficiently high to 
sample the far UV region of the extinction curve, the test 
can also be used to discriminate between different extinc- 
tion curves. In fact, for the two absorbers at Zabs '^ 2 
towards J1159-h0112 and J0016-0012 the test favours an 
extinction curve of SMC type which, at variance with the 
MW curve, yields a constant Ay from the different mea- 
surements of A(j/ — x) (bottom panels in the figure). 

In principle, the comparison of the results obtained 
from different color indices could also be used to test 
whether the absorber originates at the redshift of the 
quasar or not. We have performed this test for the cases 



in which 



-qso 



is significantly larger than Zabs, namely 



J1323-0021, J2340-0053 and J0121+0027. In all these 
cases we find that a MW-type extinction curve does not 
pass the test {Ay shows a large scatter in different color 
indices) , but an SMC extinction curve is still a viable pos- 
sibility [Ay is approximately constant). We conclude that 
this test, taken alone, is unable to rule out the possibil- 
ity that the reddening originates at the redshift of the 
quasar. This in turn implies that the reddening detection, 
even if confirmed in different color indices, must be accom- 
panied by a spectroscopic analysis of the absorption lines, 
in order to establish the location along the line of sight of 
the absorber responsible for the reddening. This conclu- 
sion is reinforced by the analysis of the quasars without 
absorption lines in our control samples: by choosing very 
reddened quasars in these samples one may obtain an ap- 
proximately constant Ay adopting an arbitrary value of 
Zahs < ^qso in the conversion from reddening to extinction. 

3.3. Extinction curve and average Ay 

In order to obtain a value of Ay for each absorption system 
of our sample we adopted one of the two types of extinc- 
tion curve (SMC or MW) and then averaged the extinction 
obtained from different color indices. In practice, we aver- 
aged the independent measurements of Ay obtained from 
two colors without common passbands. Local systematic 
errors, such as the contamination of the quasar contin- 
uum by emission lines, arc minimized in this way. In most 
cases, we averaged the measurements obtained from the 
{g — i) and (r — z) colors. For the quasars with u band 
not contaminated by Lya emission/forest (zqso < 1.67) 
we used the {u — i) and {g — z) colors, which have a larger 
spectral leverage for the detection of the reddening. For 
J2234-I-0000 (zqso = 3.015) we only had one uncontami- 
nated color index, i.e. the (r — z). 

In order to choose the extinction curve we analysed the 
SDSS spectra of the quasars of our sample in a search for 



the redshifted 217.5 nm extinction bump. Since this fea- 
ture is present in the MW curve but absent in the SMC 
curve, it can be used to discriminate between the two types 
of extinction. Within the limits of the wavelength cover- 
age of the SDSS spectra we were able to search for the 
redshifted bump in 8 quasars of our list. 

In 7 of these quasars we do not find evidence of the 
bump from the comparison of the quasar spectrum with 
the MW-type extinction curve (CCM with Ry = 3.1) 
plotted in a scale 10~5[-'*/(i+2aba)] -v^hich, apart from a con- 
stant factor, represents the transmission of the DLA sys- 
tem (see Eq. IA.6I in Appendix A) . One example of neg- 
ative detection is shown in Fig. |3i, where one can see 
that the continuum of J0016-0012 changes smoothly, as 
expected by a power law, in the region where the bump is 
expected to produce a broad, deep absorption. We cannot 
exclude that, in exceptional cases, a particular combina- 
tion of quasar emission lines may conspire to hide a dust 
bump, if present. However, it is unlikely that this hap- 
pens systematically in the 7 cases considered here, char- 
acterized by different combinations of emission and ab- 
sorption redshifts. In particular, we do not find the bump 
in J1159-h0112 and J0016-0012, the two cases for which 
the test of Fig. |21 indicates an SMC type extinction. In 
these two cases, therefore, we have two tests, one based 
on photometry and the other on spectroscopy, that con- 
sistently suggest an SMC-type extinction curve. 

Only in one case, J0121-I-0027, we found tentative 
evidence of the bump. This case is shown in Fig. Ob, 
where one can see that the continuum of J0121-I-0027 
runs smoothly at A > 600 nm, but shows an absorption 
at 400 nm < A < 600 nm. The position, width and inten- 
sity of this absorption feature are broadly consistent with 
those predicted for the bump, given the uncertainties in- 
duced by the quasar emission lines. For this absorber we 
consider the possibility that the extinction curve may be 
of Milky- Way type, without ruling out completely a SMC- 
type curve (see Table 3). We refer to Wang et al. (2004) 
for a detailed spectroscopic analysis of this case. 

For the absorbers of our list without spectral coverage 
of the bump we adopted an SMC curve. This choice is 
justified by the results of Wild & Hewett (2005) and York 
et al. (2006) which have found no evidence of the bump 
in their average absorber spectra. In any case, the results 
for an SMC and MW curve are often very similar. For the 
DLA system in front of J1323— 0021, the only reddened 
quasar of our list without adequate spectral coverage'^ of 
the bump, we expect a difference of only 0.08 dex in Ay . 

The resulting values of extinction are listed in Table 3. 
The large upper error bar that we assign to J0121+0027 
reflects the uncertainty in the choice of Ry , since a value as 
high as Ry ~ 5.5 has been proposed by Wang et (2004) for 



■^ If present, the bump of the absorber at Zabs ~ 0.716 in 
J1323-0021 should fall at the violet edge of the SDSS spec- 
trum. At shorter wavelengths, the Lya emission of the quasar 
and the Ly a forest will make difficult the search for the bump 
even in spectra with better coverage. 
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Fig. 3. SDSS spectra of two quasars with a foreground DLA system at redshift Zabs- Smooth lines: MW-type extinction 
curve {Rv — 3.1) plotted in a scale 10~f[''*/(^+^''i'=)l. See Section 3.3. The predicted positions of some emission lines 
typical of quasar spectra are indicated. Intensities are in arbitrary units. 



the absorber at Zabs — 1-388 (see footnote h). If Ry is so 
high, the extinction curve becomes flat and the conversion 
from reddening to Ay unreliable, as we explained above. 
In any case the value of extinction inferred by Wang et al. 
for this system, XogAy — +0.10 dex, lies well within our 
upper error bar. 



4. Extinction versus metal column density 

On the basis of the of the selection process of our sample 
(Section 2) we assume that the DLA systems of Table 1 
are the major source of reddening of their quasars and we 
investigate the relation between the quasar extinction and 
the DLA metal column density. The rest-frame extinction 
scales with the dust-phase column density of iron, N-pc 
[cm~^], according to the relation 



Av^ 



{41 



Nfc 



(3) 



derived in Appendix B. The term (sy^) is the mean op- 
tical cross section of the dust grains in the V band per 
atom of iron in the dust. The extinction will change be- 
tween different lines of sight tracking the variations of the 
dust-phase column density of iron and of the mean cross 
section (sy°). Our goal is to probe the behaviour of (sy ) 
in different H i regions from a simultaneous measurement 
of Ay and N-p^. The measurement of iVpe in extragalactic 
clouds is discussed in the next section. 



4.1. The dust-phase column density of iron 

The dust-phase column density of iron cannot be mea- 
sured directly from optical/UV spectroscopic data. With 
this type of observations, however, we can infer N-pc from 
the measured gas-phase column densities of zinc and iron. 
We use zinc as gas-phase tracer of iron, since zinc is 
a volatile element with little affinity to dust (Savage & 



Sembach 1996) and, at the same time, the zinc/iron abun- 
dance ratio is approximately constant, in solar proportion, 
in a large fraction of Galactic stars over a very wide range 
of metallicities (Mishenina et al. 2002, Gratton et al. 2003, 
Nissen et al. 2004). 

If zinc were completely undepleted and the intrinsic 
Zn/Fe ratio in DLA systems perfectly solar, the Zn/Fe 
ratio observed in the gas would give an exact measure of 
the iron depletion Sfc = \og{NFo/Nzn) — log (Fe/Zn)Q. In 
this case the dust-phase column density of iron would be 



iVpo = /fc Nzn 



Znjp 







(4) 



where /fc = 1 — lO''^" is the fraction of iron in the dust. 
With this expression Npc can be directly estimated from 
the measured column densities A^zn and A^fo- 

A closer inspection of interstellar data, however, in- 
dicates that zinc also can be depleted, with an extreme 
value as high as —0.6 dex in dense clouds such as that in 
front of COph (Savage & Sembach 1996). Therefore, we 
must take into account the possibility that zinc may also 
be depleted (in small amounts) in DLA systems. 

In addition, a closer analysis of stellar data shows that 
the Zn/Fe ratio tends to increase with decreasing metallic- 
ity. The effect is strong when [Fe/H] < —2 dex, attaining 
a value as high as [Zn/Fe] ~ -1-0.5 dex at [Fe/H] ~ —4 
dex (Cayrel et al. 2004). Luckily, this interval of metallici- 
ties is well below the typical metallicities of DLA systems. 
However, a mild excess has also been found in a fraction 
of Galactic stars which are less metal deficient (Prochaska 
et al. 2000, Chen et al. 2004). For stars of the thin disk the 
typical excess is [Zn/Fe] ~ +0.1 dex at [Fe/H] — —0.6 dex 
and vanishes at higher metallicities (Chen et al. 2004). 

By allowing a fraction /z„ of zinc to be in the dust 
and the intrinsic Zn/Fe ratio of the absorbers, (|^) , to 
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Table 3. Extinction and metal column densities of the DLA systems of Table 1. 



SDSS 


Zabs 


(logAv)" 


(logAv)^„d 


logiV(Znn)'= 


[Zn/Fc]" 


(logiVpo)'' 


(logiVpo)" 


(logiVFo)-^ 


J0013+0004 


2.025 


< -1.02 


-1.89, 


-1.67 


12.25 ±0.05 


0.05 ±0.07 


< 15.00' 


< 15.00" 


< 14.80" 


J0016-0012 


1.973 


-0.80+-- 


-1.22, 


-1.00 


12.82 ±0.04 


0.84 ±0.05 


15.61 ±0.04 


15.62 ±0.04 


15.50 ±0.04 


J0121+0027'' 


1.388 


-o.wtlll 


— 




> 13.32 


> 1.48 


> 16.17 


> 16.44 


> 16.19 


J0938+4128 


1.373 


< -0.69 


-1.85. 


-1.63 


12.25 ±0.05 


0.29 ±0.11 


14.80 ±0.14 


14.80 ±0.14 


14.56 ±0.14 


J0948+4323 


1.233 


< -0.51 


-1.03, 


-0.80 


13.15 ±0.01 


0.45 ±0.02 


15.82 ±0.02 


15.82 ±0.02 


15.65 ±0.02 


,11010+0003 


1.265 


< -0.88 


-0.98, 


-0.76 


13.15 ±0.06 


0.75 ±0.08 


15.93 ±0.06 


15.93 ±0.06 


15.80 ±0.06 


J1107+0048 


0.741 


< -0.58 


-1.01, 


-0.79 


13.03 ±0.05 


0.37 ±0.08 


15.65 ±0.08 


15.65 ±0.08 


15.46 ±0.11 


J1159+0112 


1.943 


-0.85l«:- 


-1.25, 


-1.02 


13.09 ±0.08 


0.45 ±0.11 


15.76 ±0.10 


15.76 ±0.10 


15.59 ±0.12 


J1232-0224 


0.395 


< -0.50 


-1.10, 


-0.87 


12.93 ±0.12 


0.83 ±0.16*^ 


15.72 ±0.12 


15.72 ±0.12 


15.60 ±0.13 


J1323-0021 


0.716 


-0.36tr,l 


-0.56, 


-0.33 


13.43 ± 0.05 


1.14 ±0.06 


16.26 ±0.05 


16.28 ±0.05 


16.16 ±0.05 


J1501+0019 


1.483 


< -0.79 


-0.92, 


-0.70 


13.10 ±0.05 


0.40 ±0.09'' 


15.74 ±0.08 


15.74 ±0.08 


15.56 ±0.11 


J2234+0000 


2.066 


< -0.61 


-1.60. 


-1.38 


12.46 ±0.02 


0.57 ±0.06 


15.19 ±0.03 


15.19 ±0.03 


15.05 ±0.04 


J2340-0053 


1.360 


"■°'-0.32 


— 




12.62 ±0.05 


0.54 ±0.13 


15.33 ±0.07 


15.33 ±0.07 


15.18 ±0.10 



" Measured from the observed reddening; using Eg. jA.lOl for the convcrsJQn to rest frame extinction. 

Extinction predicted from Eq. (11) of Vladilo & Pcroux (2005) for G — 0.6 and G — 1, respectively. 
^ See references in Table 1. 

Calculated assuming Zn completely undepleted and [Zn/Fe]a, — 0. 
'^ Calculated using a scaling law of interstellar depletions (Vladilo 2002b; ezn — 1^ Vzn — 42.5) and [Zn/Fe]a — 0. 
■'^ As in the previous column, but with [Zn/Fe]a — +0.1 dex. 
^ Dust-free system; upper limit computed introducing logA''(Znll) < 12.40 and log A''(Fell) > 14.91 in Eqs. (^ and 1^ . 

Central and upper values of (logAv) calculated using a CCM model with R^ — 3.1 and R^ — 5.5, respectively. Lower value calculated using 
the SMC extinction curve (see Section 3.1). 

A'^(Ni) or A'^(Cr) used as a proxy of A'^(Fe), which is not measured; Ni and Cr are good tracers of Fe in stars with [Fe/H] > — 2 dex (e.g. Ryan 
et al. 1996); Ni, Cr and Fe have very similar interstellar depletions (e.g. Vladilo 2002a); we adopt [Ni/Fe]— and [Cr/Fc]— in the conversion. 



deviate from the solar value, we obtain the more general 
expression 



A^i 



Nv, 



Fe 



/Fe 



(1 - /zn) 



Zn 
Fe 



(5) 



With an educated guess of (Zn/Fe)a it is possible to es- 
timate A'fc from the column densities iVzn and A'fc also 
in this general case, assuming that fzn scales with /pe 
according to a scaling law of interstellar depletions cali- 
brated in the ISM. Details on this method can be found 
in Vladilo (2002a, 2002b). 

In Table 3 we give the dust-phase column densities 
of iron calculated for three different cases: (1) zinc com- 
pletely undepleted and [Zn/Fe]a — (column 7); (2) zinc 
depletion scales with iron depletion and [Zn/Fe]a = (col- 
umn 8); (3) zinc depletion scales with iron depletion and 
[Zn/Fe]a = -HO.l (column 9). 

The comparison between columns 7 and 8 indicates 
that the results are almost unaffected by the possible de- 
pletion of zinc. Only when the depletion level is particu- 
larly high, i.e. when [Zn/Fe] > 1 dex, the correct relation 

(O yields slightly higher values of iVpo than the approxi- 
mate relation Q). 

The comparison between columns 8 and 9 indicates 
that a mild excess of the intrinsic abundance ratio 
[Zn/Fe] a yields a similar mild decrease of Npc with re- 
spect to the case [Zn/Fe]a — 0. The higher the observed 
[Zn/Fe] ratio, the lower the difference. For the quasars 
with reddening detection the observed [Zn/Fe] ratio of the 
DLA happens to be relatively high (column 4 in Table 3). 
Therefore the results of the present work are weakly af- 



fected by an excess of [Zn/Fe] a of the type observed in 
stars of moderately low metallicity. 

In fact, for two of the quasars with reddening detec- 
tion (J0121-I-0027 and J1323-0021) the metallicity of the 
intervening DLA system is relatively high (Wang et al. 
2004, Peroux et al. 2006) and we do not expect an excess 
of [Zn/Fe] a given the trend of [Zn/Fe] observed in Galactic 
stars (Chen et al. 2004). The absorber with lowest metal- 
licity among the 5 quasars with reddening detection is 
that at Zabs = 1.944 towards J1159+0012 ([Zn/H] ~ -1.4 
dex). This is the absorber in which [Zn/Fe]a could be, in 
principle, most enhanced. As we show below, the results 
of this paper are unaffected even in this case. 

4.2. The empirical relation between Ay and A^pe 

In Fig. 0] we plot the rest frame extinction. Ay, versus the 
dust-phase column density of iron, N-pc, for different types 
of absorbers of the local and high redshift Universe. 

The data obtained from our reddening detections at 
the 2 a level are indicated with circles. The filled circles 
are the quasars in which the DLA system of Table 1 is the 
major source of reddening. The empty circle is J2340-0053, 
for which an additional contribution to the reddening may 
be present. For J0121-I-0027 we also plot the extinction 
derived by Wang et al. (2004) from the comparison with 
composite SDSS quasar spectra (filled square). 

For the cases without reddening detection we plot a 
2 a upper limit. The values of N-pc are estimated from Eq. 
(0), case [Zn/Fe]a=0. The implications of these choices are 
discussed below. 

The triangles in the figure represent two quasars with 
photometric reddenings and DLA metal column densities 
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Fig. 4. Extinction, Ay, versus dust-phase column density of iron, Nfb, in QSO absorbers and interstellar clouds. 
The extinction is calculated in the rest frame of the absorber. Circles and arrows pointing down: Ay measurements 
and upper limits from this work; empty circle: J2340— 0053. Triangles: QSO/DLAs from Ellison et al. (2005); empty 
triangle: BG438— 436. Square: extinction for JG121+0027 from Wang et al. (2004). Crossed diamonds: sub-samples of 
Mgii absorbers from York et al. (2006). Open diamonds: sub-samples of Can absorbers from Wild et al. (2006). Open 
star: SMC line of sight towards Sk 155 (Welty et al. 2001, Fitzpatrick 1985). Dots and crosses: Milky Way interstellar 
lines of sight from Jenkins et al. (1986) and Snow et al. (2002), respectively. Strip: ztlcr interval around the linear 
regression of the MW data. Dashed line: linear regression of the MW data with fixed unit slope. See Sect. 4.2. 



obtained in the framework of the CORALS survey (Ellison 
et al. 2005). These quasars are not included in SDSS cat- 
alog. To estimate Ay we adopted the normalized colour 
{B—K)n in Table 3 of Ellison et al. and an SMC extinction 
curve. Only two quasars of this survey (B0438— 436 and 
B0458— 020) have a significant {B — K) reddening and, at 
the same time, Zn and Fe column densities (Akerman et 
al. 2005). The quasar B0438— 436 is represented with an 
empty triangle, rather than a filled one, because its red- 
dening seems to be contaminated by dust in the quasar 
environment (Ellison et al. 2005). 

The crossed diamonds in the figure are taken from a 
study of the average reddening properties of 809 Mg ii sys- 
tems in SDSS quasar spectra (York et al. 2006). The red- 
shift interval covered by this study is similar to that of the 
DLA systems of Table 1. However, Mgii absorbers may 
have, in general, H i column densities well below the DLA 
regime. For comparison with our data we have consid- 
ered the sub-samples of York et al. with largest equivalent 



widths of Mg II because in these sub-samples the fraction 
of DLA systems is expected to be the highest. The re- 
sults shown in the figure represent the mean values of Ay 
and A'fc for the sub-samples S8 (251 systems with equiv- 
alent widths > 2 A) and S26 (97 systems with equivalent 
widths > 2.5 A and A{g — i) < 0.2 magnitudes). Sub- 
samples with saturated Zn ii and Fe ii mean profiles were 
not considered. The mean extinction was estimated from 
the mean E{B — V) of these sub-samples using the SMC 
value Ry ~ 2.74 (Gordon et al. 2003). The mean A^fc was 
derived by inserting the mean A^(Znii) and [Zn/Fe] of the 
sub-samples in relation @. 

The open diamonds are taken from the study of Can 
systems at z --^ 1 in SDSS quasar spectra by Wild et 
al. (2006). In order of increasing extinction, they repre- 
sent the mean values of Ay and A'fc for the sub-samples 
labelled Low-Wa3935, All, and High-WAsgas, respectively. 
We estimated these mean values from the corresponding 
mean values of E{B — V), N{Znu) and [Zn/Fe] given by 
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Wild et al. (2006) for the 27 absorbers analysed for el- 
ement column densities. Also in this case we assume an 
SMC type of dust. 

The quasar absorption data in Fig. 21 are consistent 
with an increase of Ay with increasing Np^. A linear trend 
between these two quantities is predicted by relation Q 
if the mean extinction per atom of iron is approximately 
constant in Hi interstellar regions. However, the number 
of data points for the individual quasar absorbers (circles 
and triangles) is too little for deriving a statistical corre- 
lation. In order to cast light on the behaviour of high red- 
shift clouds we added in the same figure results obtained 
for local interstellar clouds. 

The dots and crosses in the figure represent Milky- 
Way lines of sight with available measurements of Ay, 
7V(Htot) = A^(Hi) -t- 2iV(H2) and iV(Feii), taken from 
Jenkins et al. (1986) and Snow et al. (2002), respectively. 
For these data we derived iVpc assuming that the local in- 
terstellar abundance of iron is solar. The error budget in 
the estimate of Npc is dominated by the error of iV(Htot), 
which is typically of ^ 0.1 dex. The extinctions were taken 
from the catalog of Neckel & Klare (1980), under the en- 
try Av{MK). Uncertainties of these Ay values are in the 
order of 0.1 mag. Values with Ay < 0.1 mag were not 
considered. Lines of sight with iV(Hi) below the defini- 
tion threshold of DLA systems and with saturated iron 
lines were rejected. 

The Milky- Way data are highly correlated, with linear 
correlation coefficient r = 0.95 and slope m = 0.81 ± 0.27. 
The strip in the figure represents the dispersion of the 
data around the linear regression. The dashed line is the 
linear regression with fixed m = 1, a slope consistent with 
the free regression analysis. 

The star symbol in the figure represents the SMC star 
Sk 155, the only Magellanic sightline for which we could 
determine Ay and Npc- The extinction was estimated us- 
ing the SMC value of E{B-V) from Fitzpatrick (1985) 
and the mean SMC value Ry = 2.74 (Gordon et al. 2003). 
We estimated A^fo from relation (jSJ using the Zn and Fe 
column densities of Welty et al. (2001) integrated over 
the SMC radial velocities. In spite of the differences in 
metallicity level and type of extinction curve relative to 
the Milky Way, the data point of Sk 155 is in excellent 
agreement with the Milky- Way data. 

Most of the DLA measurements (filled circles and tri- 
angles) follow the trend of Milky- Way interstellar data. 
Uncertainties in the choice of the extinction curve are un- 
likely to alter this result. Even in the case of J0121-I-0027, 
for which these uncertainties are large, the DLA data point 
seems to be consistent with the interstellar trend (the 
value of -/Vpc is a lower limit owing to saturation of the 
Znii lines). 

Also the uncertainty in the adopted value of [Zn/Fe]^ 
in Eq. Q is unlikely to affect the agreement between 
DLA systems and interstellar clouds. As we said above, 
J1159-I-0112 is the most likely candidate for a mild en- 
hancement of [Zn/Fc]^ among our quasars with reddening 



detection. By adopting an enhanced [Zn/Fe]a in this case 
(last column of Table 3), we still find a good agreement 
with the Milky- Way data. 

In spite of this general agreement, an excess of ex- 
tinction relative to the Milky- Way trend is found in 
J2340-0053 (empty circle) and B0438-436 (empty trian- 
gle). Quite interestingly, these are exactly the cases for 
which we suspect the presence of extinction sources in ad- 
dition to the DLA system. Therefore, the excess may be 
due to the presence of additional absorbers rather than to 
a deviation of individual DLAs from the interstellar trend. 

Also the Ay upper limits that we derive are gener- 
ally consistent with the interstellar data. This conclusion 
would be strengthened by adopting an enhanced [Zn/Fe]o 
in Eq. ©, instead of [Zn/Fe]^ = 0, since in this case the 
limits would shift slightly to the left (see Table 3). 

In summary, the data collected in the figure suggest 
that the trend between Ay and N-p^ is remarkably similar 
in interstellar clouds of the Milky Way and the SMC and in 
DLA systems with different metallicities (from ~ — 1.4 dex 
up to ci -1-0.5 dex relative to solar) and redshifts {z ~ 0.7 
up to z ~ 2). Also Can absorbers and the Mgii absorbers 
with highest values of equivalent width follow the same 
trend. 

In the interpretation of the results of Fig. ^ one must 
take into account that some parts of the plot are not ac- 
cessible to the observations. For instance, if Ay is too 
large, the background quasar drops out of the SDSS sam- 
ple. If it is too low, the reddening cannot be detected. 
The detection limit, in addition, will vary as a function 
of the redshift and of the scatter of the intrinsic colors 
of the quasars. These effects may alter the distribution of 
detected DLA systems in Fig. 0] and, in principle, could 
induce an artificial trend among the data points. 

The agreement between the different sets of high red- 
shift data in Fig.^lis impressive, considering the different 
methods employed for measuring the reddening. Also the 
agreement between high redshift and local data is remark- 
able, considering the different methods used to derive Ay 
and A'fc (for instance, the educated guess of [Zn/Fe]a was 
applied to derive A^fc at high redshift, but not in the Milky 
Way) . The general agreement between DLA and interstel- 
lar data speaks against the existence of an artificial trend 
induced by selection effects of DLA detections. The impli- 
cations of this general agreement are discussed below. 

5. Discussion 

5.1. Implications for the properties of dust grains 

The existence of a common trend between Ay and Npe 
in Hi regions of galaxies with different metallicities and 
redshifts indicates that the dust grain parameter 



(6) 



(.-) = 1.007x10-- ^l^^^^i^^l^ 

representing the mean extinction in the V band per atom 
of iron in the dust (see Appendix B), is approximately 
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constant in a large variety of neutral regions. The scatter 
of the Milky Way data is of only 0.16 dex at 1 cr level. The 
few measurements in DLA systems do not show, so far, 
evidence for a larger dispersion. 

In principle, the V band could be less sensitive than 
other spectral bands to variations of the extinction prop- 
erties of the grains. Notwithstanding, the small scatter of 
(s^) is rather remarkable given the fact that the value 
of this parameter is determined by at least 4 different 
properties of the grains: their size, a, internal density, 
g, abundance by mass of iron inside the grain, Xpe, and 
extinction efficiency factor, Q\ , which is related to the 
geometrical and optical properties of the grains (Spitzer 
1978). We know that these parameters vary among inter- 
stellar clouds. In particular, the variations of the extinc- 
tion curves are attributed, in large part, to variations of 
the grain size distribution (Draine 2003). Therefore, in or- 
der to produce an approximately constant (s^) there must 
be a physical mechanism able to compensate the variations 
of the individual parameters that appear in the above ex- 
pression. We tentatively propose the following mechanism 
based on heuristic considerations. 

We consider two types of interstellar environments: (1) 
regions where grain destruction mechanisms are efficient 
and (2) clouds where the grains are protected from de- 
struction. In the second case we expect that the grain 
size a can be larger, on the average, and the volatile el- 
ements more easily incorporated into the grains than in 
the first case. The relative abundance of iron inside the 
grains, Xpo, must be lower in the second case, when also 
volatile elements are incorporated in the grains, than in 
the first case, when refractory elements, such as iron, are 
dominant. Since the most abundant volatile elements (e.g. 
carbon) are lighter than the most abundant refractory ele- 
ments (e.g. iron) we expect, in addition, that in the second 
case the mean density g is lower than in the first case. 

Combining all together, we expect the following vari- 
ations of the grain parameters to occur passing from a 
grain-destructive to a non-destructive environment: an in- 
crease of a and, at the same time, a decrease of g and Xpc. 
From relation © we expect that (sy°) ~ (a gXpc)~^ since 
(Tg w a^ and Vg ^ a"^. Given this dependence of (s^), the 
changes of a could compensate the changes of g and Xpc- 

A much larger database of interstellar data at low and 
high redshift are required to understand the viability of 
the simple mechanism of compensation proposed here. 
Contraints on this mechanism could be obtained by study- 
ing the slope of the Ay and A^pe relation and analysing 
possible variations of (sy°) versus Ry- Given the paucity 
of quasar absorption extinction measurements these stud- 
ies are still premature at high redshift, but could be the 
subject of future investigations. 



5.2. Implications for the dust obscuration bias 

The extinction of the quasars due to the intervening DLA 
systems is expected to produce a selection effect in opti- 



cal, magnitude-limited surveys: the absorbers with highest 
values of extinction would be systematically missed from 
the surveys as a consequence of quasar obscuration (Fall 
& Pei 1989). 

Evidence for this effect has been searched by com- 
paring the statistics of optical surveys and radio-selected 
surveys, the latter being unaffected by the extinction. A 
marginal signature of the bias has been found in this way 
(Ellison et al. 2001, Akerman et al. 2005), but the statis- 
tics are still based on small numbers and the estimate of 
the effect uncertain due to the difficulty of observing all 
the quasars of the radio-selected sample, independently of 
their magnitude, in the optical follow-up. 

In a previous work we have proposed a new ap- 
proach for quantifying the effect of quasar obscura- 
tion on the statistics of DLA systems. Using a rela- 
tion between the extinction and the metal column den- 
sities of the absorbers we invert the frequency distri- 
butions of H I column densities and metallicities mea- 
sured from magnitude-limited surveys and derive the un- 
biased frequency distributions (Vladilo & Peroux 2005). 
The existence of a general relation between extinction 
and metal column density is fundamental in this type 
of approach. The results of the present work lend quan- 
titative support to the relation adopted in our previ- 
ous work^, Av ~ 1.85 x 10"" G/pc(Z) A^h ^, which was 
mostly based on interstellar data. In Table 3, Col. 4, we 
list the values of extinction predicted from this previous 
relation for the DLA systems of the present work. One can 
see that the predictions are in most cases consistent with 
the measurements or upper limits. The agreement with 
the measurements is better for G = 1 than for G — 0.6 (in 
this latter case the extinction is mildly underestimated). 
These results indicate that the quantitative predictions of 
the obscuration bias presented by Vladilo & Peroux (2005) 
are realistic and, in any case, not overestimated. 

In order to confirm the general validity of the relation 
between extinction and metal column densities in DLA 
systems, it would be important to obtain more measure- 
ments at higher redshifts since the present sample is lim- 
ited at Zabs — 2, while a large fraction of DLA systems is 
observed to be at higher redshift. In particular, it will be 
crucial to calibrate the relation at high values of extinction 
{Ay > 1 mag), when the obscuration bias is most critical. 
The slope m ~ 0.8 of the interstellar data is consistent 
with a mild decrease of the extinction per atom of iron 
in the dust with increasing metal column density (strip 
in Fig. 0}. For an exact estimate of the obscuration bias 
it is critical to verify the existence of this effect in DLA 
systems. 



" See Eq. (11) in Vladilo & Peroux (2005), where 
Z = Nzn/NH and G = 1 (0.6) for MW (SMC) type 
of dust; £,{Xv) = 1 in the rest frame of the absorber; 
fFc{Z) = i {arctan ( * " q 4 ' "* ) -I- §} . Apart from a normal- 
ization factor, the parameter G is equivalent to {s^) for the 
special case of a single family of identical, spherical grains. 
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Fig. 5. (a) Extinction, ^y, versus total hydrogen column density, iVn, of DLA systems and interstellar clouds. Same 
symbols as in Fig.|3| Dashed line: mean MW Ay/N-a ratio for Ry = 3.1 (Bohhn et al. 1978). Dotted line: mean SMC 
A-v /N-n ratio (Gordon 2003). (b) Dust-to-gas ratio, Ay/Nn, versus metallicity for the same set of data. The labels 
indicate the regions of the diagram occupied by MW and SMC interstellar clouds. 



5.3. The dust-to gas ratio Ay/Nn 

In Fig. E^ we plot Ay versus the total hydrogen column 
density for the same data set of Fig.01 For DLA absorbers 
we assume that Nn — A^(Hi), i.e. we neglect the molec- 
ular contribution. Two of our systems with reddening de- 
tections are not plotted because they do not have a direct 
determination of iV(H i) (see Table 1). This is also the case 
for the Mgii absorbers of York et al. (2006) and the Can 
absorbers of Wild et al. (2006) discussed in Section 4.2. 

The comparison with the interstellar data shows that 
DLA systems do not follow the typical Milky Way relation 
of Bohlin et al. (1978), shown as a dashed line in the figure. 
Most systems lie below the MW trend and the addition of 
a molecular contribution to the hydrogen column density 
would not change this conclusion. 

The analysis of the data versus metallicity reveals the 
existence of a regular behaviour. This can be seen in Fig. 
Oj, where we plot the dust-to-gas ratio Ay/iVn versus 
metallicity [Zn/H] for the same systems. In spite of the 
limited size of the DLA sample, the data suggest the ex- 
istence of a trend of increasing dust-to-gas ratio with in- 
creasing metallicity. The trend nicely fits the regions of the 
diagram occupied by MW and SMC interstellar clouds, 
suggesting its general validity at low and high redshift. 

The observed trend can be interpreted as a natural 
consequence of the approximate constancy of (s^) dis- 
cussed above. In fact, dividing both sides of Eq. JSJ) by 
TVh and inserting relation Q, we obtain 






/fc A^zn /ZnV 
(l-/z„) A^H UeA 



(7) 



If (sy) is constant, we expect from this expression a linear 
trend between dust-to-gas ratio and metallicity (Zn/H), 
modulated by variations of the iron and zinc depletions; 
the intrinsic abundance ratio (|^) is not expected to 
show strong variations. 

If the trends of Fig. 0] and |5l will be confirmed by a 
larger set of data, then it will be natural to explain the 
vahdity of the Bohhn et al. (1978) MW relation, Av/Nk - 
constant, as the consequence of the constant metallicity 
level of MW clouds in conjunction with a constant (sy ). 

The deviations of the DLA dust-to-gas ratios from the 
Bohlin et al.'s value indicate that the MW ratio should 
not be applied to obtain indirect estimates of extinc- 
tions or A^H in extragalactic research. By taking into ac- 
count the metallicity and the extinction curve of the ab- 
sorbers one can, however, obtain reasonable estimates. 
Since quasar absorbers are, in general, metal deficient, the 
typical Av/Nn ratio of the SMC (dotted hne in Fig.EJ is 
more appropriate than the MW ratio for indirect estimates 
of extinctions or A^h- The choice of an SMC dust-to-gas 
ratio is in line with the lack of 217.5 nm bump in most 
quasar absorbers (York et al. 2006). 

6. Summary and conclusions 

We have implemented a technique for measuring the red- 
dening of quasars using the photometric and spectroscopic 
SDSS database. For each quasar under investigation we 
build up a control sample to analyse the distribution of the 
colors at a given redshift and magnitude. Quasars show- 
ing absorption lines redwards of the Ly a emission in their 
spectra are rejected from the control sample. The median 
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and the dispersion of the color distribution are used to 
estimate the reddening and its uncertainty. To minimize 
systematic errors related to the presence of quasar emis- 
sions in particular passbands, the reddening is measured 
in different color indices. 

We have applied this technique to the complete sam- 
ple of 13 SDSS quasars for which a single intervening DLA 
system with zinc absorption lines had been previously de- 
tected in spectra of high (or intermediate) resolution. Most 
of the spectra of these quasars do not show evidence of 
potential sources of reddening other than the DLA sys- 
tems. We detect reddening at the 2 a level in five of these 
quasars. In each case the detection is confirmed in different 
color indices. The comparison of the color excess measured 
in different pairs of bandpasses is consistent with an origin 
of the reddening in the DLA systems. To our knowledge 
this is the most direct evidence of DLA extinction up to 
z ~ 2 obtained so far. 

In order to discriminate the extinction curve of the ab- 
sorbers we have compared the reddening in different color 
indices and investigated the region of the SDSS spectra 
where the redshifted 217.5 nm MW-extinction bump is 
expected to fall, if present. An SMC-type extinction curve 
is generally consistent with the available data, in agree- 
ment with previous studies of the extinction curves of Ca ii 
and Mgii quasar absorbers (Wild & Hewett 2005, York et 
al. 2006). The only possible exception is J0121+0027, for 
which we cannot exclude that the MW bump is present, 
as claimed by Wang et al. (2004). 

After converting the quasar reddening into V^-band ex- 
tinction in the rest frame of the DLA system. Ay (see 
Appendix A), we have investigated the relation between 
Ay and the dust-phase column density of iron of the ab- 
sorber, N-pc- Our measurements and upper limits are con- 
sistent with a rise of Ay with increasing A'po, but the 
sample is too small to perform a correlation study. 

By comparing these data with measurements of Ay 
and TVpe m Milky Way and SMC interstellar clouds, we 
obtain the main result of our work: the high-redshift data 
follow remarkably well the trend of the local interstellar 
data, which show a linear correlation log Ay - logA^Fo- 
Only in one quasar, J2340— 0053, do we find an excess 
of extinction, possibly due to neutral gas at a redshift 
different from that of the DLA system. 

Consistent results are found from the analysis of the 
two quasars with {B — K) reddening detections recently 
reported by Ellison et al. (2005): one DLA system with 
Znii lines follows the Milky- Way trend and another one, 
towards B0438— 436, shows an excess of extinction. In this 
case the excess is likely to originate in dust of the quasar 
environment (Ellison et al. 2005). 

Finally, from the studies of Wild et al. (2006) and York 
et al. (2006), we find evidence that also the Can quasar 
absorbers at z ~ 1 and the Mgii absorbers with highest 
equivalent widths, statistically more similar to DLA sys- 
tems, follow the interstellar trend between extinction and 
dust-phase metal column density. 



The existence of a linear relation between Ay and 
A^Fo shared by high redshift DLA galaxies and by local 
clouds of the Milky Way and the SMC suggests that the 
mean extinction per atom of iron in the dust, (sy°), is 
approximately constant in galaxies with different levels of 
metallicity (from [Zn/II]~ —1.4 dcx up to ~ -1-0.5 dex) 
and look-back times (up to ^ 10 billon years before the 
present). This result in turn suggests the existence of a 
mechanism which tends to compensate variations of (sy?) 
resulting from changes of the dust grain properties in dif- 
ferent interstellar environments. We propose that when 
the grain size a increases, the density of the grain, p, and 
the abundance by mass of iron in the dust, Xpe, decrease 
(and vice versa). These variations are expected to occur 
when passing from a dust-destructive to a dust-protected 
environment (and vice versa). 

The existence of a well defined trend between metal 
column densities and extinction in DLA systems at dif- 
ferent redshifts has important implications for quantify- 
ing the selection effect of quasar obscuration (Fall & Pei 
1989). The results of the present analysis lend quantita- 
tive support to the estimates of DLA extinction performed 
in our previous study of the obscuration effect (Vladilo & 
Peroux 2005). 

Finally, from the analysis of the dust-to-gas ratio 
Ay /Nyi in the DLAs of our sample, we find significant de- 
viations from the Milky- Way ratio Ay /Nyi ~ 5.3 x 10~^^ 
mag cm^ of Bohlin ct al. (1978). The dust-to-gas ratio 
appears to increase with the metallicity of the absorber, 
following a trend that fits the SMC and MW data points, 
but more data are necessary to confirm this result. We 
argue that the constancy of the dust-to-gas ratio in lo- 
cal interstellar clouds is the result of their constant level 
of metallicity in conjunction with the lack of significant 
variations of {s^). The MW dust-to-gas ratio should not 
be used in extragalactic studies for indirect estimates of 
extinctions or Hi column densities. 

An important outcome of the present investigation 
stems from the comparison of the control samples be- 
fore and after the rejection of quasars with absorptions 
lines redwards of the Lya emission. After the rejection, 
the color distributions of the control samples experience a 
small but systematic blue shift indicating that these ab- 
sorption lines do contribute to the reddening of quasars. 
This result is in line with the recent detection of redden- 
ing due to Can (Wild & Hewett 2005, Wild et al. 2006) 
and Mgn (Khare et al. 2005; York et al. 2006) absorption 
systems. Even in the absence of a direct measurement of 
the H I column densities of the Ca n and Mg n absorbers, 
the results of these previous investigations provide indi- 
rect evidence that DLA systems should give a signal of 
reddening. The present investigation is in line with this ex- 
pectation, even if the sample is too small to draw general 
conclusions on the statistical properties of the extinction 
in DLA systems. 

At the present state of the observations it is hard to 
understand whether or not there is a discrepancy between 
the reddening detection in Can, Mgn, and DLA absorbers 
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at z < 2, on the one hand, and the upper hmit found by 
Murphy & Liske (2004) in DLA systems at z ~ 2.8, on 
the other. Future studies of extinction and metal column 
densities in DLA systems should be aimed at obtaining 
measurements at higher redshifts (z > 2), at high values 
of extinctions {Ay > 1 mag) and at different values of Ry ■ 
With such measurements one should be able to understand 
if the relation between extinction and dust-phase metal 
column density is indeed universal, to probe the redshift 
evolution of the dust in cross-section selected galaxies, to 
quantify with accuracy the effect of quasar obscuration, 
and to probe the mechanism that makes (s^°) relatively 
constant in interstellar clouds. 
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Appendix A: Conversion of the quasar reddening 
into rest frame extinction 

Our goal is to estimate the rest frame extinction, Ay, from 
the measurement of the quasar reddening in the observer 
frame, A(y — x). We start from the definition of the ap- 
parent magnitude of a quasar in a photometric band x, 



We now express the transmission of medium in terms 
of the extinction in the absorber rest frame, A\ (magni- 
tudes). By introducing the normalized extinction curve of 
the absorber, ^(A) = Ax/ Ay, we obtain 



5, (A) = 10 



-OAAviW 



(A.6) 



The reddening of the quasar in the color index (y — x) 
in the observer's frame is A(y — x) = A' — A'^. Combining 
the above expressions we derive 



2.5 log 



F{X)S^{X)dX 



(A.l) 



A{y-x) = 2.5 log 



J^^2 iQ-OAAv «(Ao)i(A^) S^{X) dX 



/;-L(Ae)5,(A)dA 



A.., 



where m^o is the zero point of the photometric system; 
F(X) is the monochromatic flux of the quasar received out 
of the Earth's atmosphere (erg s~^ cm"^ A~^); Sx{X) is 
the response curve of the passband x (we include in this 
term the transmission of the Earth's atmosphere and the 
response curve of the detector); Xxi and Xx2 are the lower 
and upper cutoff of the bandpass 



-2.5 log 



j^y2 io-o.iAva^a)^x^)SyiX)dX 
ltL{Xe)Sy{X)dX ■ 
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The reddening is independent of the zero point of the 
photometric system and of the cosmological parameters. 
The dependence of A(y — x) on the redshift and the spec- 



Ze IS 



F{X) 



tral distribution of the quasar is modest since the quasar 
The flux received at Earth from a quasar at redshift continuum varies smoothly in the bandpasses [A,,, A, J 

and [A.yi , Xy^]. The spectral energy distribution can be ap- 
proximated with a power law L{X) oc A"^ with spectral 
index a\ = —{2 + a^). 

To solve the above equation in Ay we perform the 
change of variables 



L{X,)Hl{l 



4^c2Z2(Ze) 



£(Ae)(l+Ze)' 



where the luminosity distance Dl = cZq{ze) (1 - 
is defined through the function (Weinberg 1972) 



(A.2) 

Ze)/Ho 



t(Aa) = -0.4 1n(10)e(Aa): 



Zq{Zi,) 



q^l + Ze) 



qoZe 



(Qo 



(l + 2qo2e)5 -l]} 



and 



Here L{X) is the monochromatic luminosity emitted by 
the quasar (erg s~^ A~^), and Ae = A/(l -I- z^). 



From this we have 
m^ = Cxize) ~ 2.5 log 



^xi 



L{Xe)SxiX)dX 



(A.3) 



r = -0.4 1n(10) A{y - x) . 
Then ea. HA.7p becomes 

Jx:^e^-'<^^^)LiX,)Sy{X)dX 
JxrLiX,)SyiX)dX 



J^^-e^-'(^-)L{Xe)Sx{X)dX 



where the term C'x{ze) depends on the zero-point of the 
photometric system, the quasar redshift and the cosmo- 
logical parameters. 

We now consider the extinction due to an intervening 
interstellar medium i that lies at the absorption redshift 
Zabs in the direction of the quasar. We call Si{X) the trans- 
mission of the medium as a function of wavelength in the 
absorber rest frame [S'i(A) = 1 for a transparent medium]. 
The reddened magnitude of the quasar is 



m; = C,(ze)-2.51og/ L{Xe)S,{Xa)Sx{X)dX, (A.4) 

where Aa = A/(l + Zabs)- 

By definition, the quasar extinction in the observer's frame 
is A' = m' — TTix ■ From this we derive the relation 



/^^^^L(Ae)5,(A)dA 



(A. 



We now express the integrals in the form of sums, cal- 
culating the functions in n^ equispaced points Xi in the 
interval [A2;j,A2;2] and Uy equispaced points Xj in the in- 
terval [Xy^^Xy^]. We obtain 
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where Si = Sx{Xi), Sj 

U = L[X^/{l + Ze)]. 

By introducing the weights 



Sy{Xj), U = t[Xi/{l + Zabs)], 
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where the term between angle brackets is an average cal- 
_i culated by integrating Q^^' Ug over the distribution of 

grain sizes a (see e.g. Pei 1992). The total extinction is 

w,{QxOg)i (B.2) 



E^-.e^'** 



(A.IO) 



A> 



1.086 A/'d E' 



where YTi=\'^x, = 1 and X]J=i "^v, = 1- 



Ki) 



The extinction Ay can be calculated from this expression where the weight Wj = M^^' jMd is the fraction of dust 

or estimated from its simplified version. An approximate grams m each family. 

value of Ay can be derived in the following way. We now use iron as a tracer of the metals in the dust. 

We call Vg and q the volume and the internal mass density 

Let m^ be the value of i for which the weight w^^ is max- of a grain, respectively, and Xpc the abundance by mass of 

imum, i.e. w^^ = max{w^.}"j;;;^ and ruy the value of j iron in the grain. The dust-phase column density of iron 

w„, is maximum, i.e. w„,,., = max|w„,)"_^i . We atoms in the j-th family of grains is 



for which Wy^ is maximum, i.e. Wr,i^ = max{wj^^.}"^P We atoms in the 



take the terms Wm^ e^vt^^ o^j^^j yj^ ^Avt^y ^^^ qJ ^Yie 
sums and obtain 
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If ti varies smoothly in the intervals [A^r^jA^^] and 
[Ayi,Ay2], so that Av \U-tmJ < 1 and Ay \tj-tmy\ < 1, 
we can expand each exponent inside the sums into a series 
e^ ~ 1 -f a;. After some algebra, we obtain 

e'- = [l + Av{{t)y-t^J] e^-'-^y x 

x{[l + Av((t),-t,„J] e■^-*'-}-^ (A.12) 

where {t).^, = YJ^l^ w^^U and {t)y = Y.'j=i '^vj^i ■ 

If \{t)x - UnA < 1 and \{t)y - trny\ < 1, then we obtain 

r ~ Avytjyiy ~ tm^) ) 
which is equivalent to 
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For narrow bandpasses (|A.13|I gives the correct solution. 
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where Apc is the atomic mass of iron. The average between 
angle brackets is calculated by integrating Vg g'^^^ X^^ 
over the distribution of grain sizes. The total dust-phase 
column density of iron is 
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Apo mn 
Combining ljB.21) and (|B.4|I we have 
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where 

(sf) = 1.007x10^22 



E,-"^j(Qa (Tg)j 



E, 



! Vg g Xpc )j 



(B.4) 



(B.5) 



(B.6) 



is the mean optical cross section (magcm^) per atom of 
iron in the dust. 



Appendix B: The mean extinction per atom of 
iron in the dust 

We consider a line of sight intersecting dust embedded in 
a neutral region. We call A^ the number of dust grains per 
cm^ along the line of sight, ag the geometrical cross section 
of a single grain and Q\ the extinction efficiency factor, i.e. 
the ratio between optical and geometrical cross section of 
the grain (Spitzer 1978). The dust will consist of different 
families of grains, each one with its own physical/chemical 
properties and, in particular, distribution of grain sizes a. 
The extinction due to the j-ili family is 



A^^ = imQNlf'' {Qxag), 



(B.l) 



